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In a recent Letter ||l| , Stallmach et al. reported pulsed 
field gradient (PFG) NMR measurements of the time- 
dependent diffusion constant D{t) in packings of water 
saturated sands. According to theory D{t) decreases 
from the bulk water value Dq with increasing observation 
time t due to restrictions imposed by the pore surface: 
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Here S/V is the surface-to- volume ratio of the pore space 
and r — ^/D^t is the diffusion length. Stallmach et al. 
studied samples with different grain diameters dg and 
found that S/V oc d~^-'' , which they interpreted with a 
fractal picture. If dg is identified as the upper cut-off 
scale L for a fractal surface of dimension , one expects 
S/V cc d^^^^. In this Comment, we argue that the 
analysis of Ref. ^is flawed and we propose an alternative. 

The key point is that Eq. ^ was derived for a nonfractal 
surface such that S/V is constant as the length scale 
r varies. In Ref. Eq. |l| is used for 2 fim < r < 
10 /im while dg is in the range 100-1000 /im. For the 
analysis to be valid, the surface would have to be smooth 
below 10 /im and abruptly turn into a fractal above this 
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scale. This is implausible as sands and rocks are known 
to have fractal surfaces below 1 fim |^. For example, 
BET measurements on rocks (for which the measurement 
scale is r R:! 0.4 nm) yield S/V values one to two orders 
of magnitude greater than PFG NMR results [|| . 

To show how Eq. ^ is modified for fractal sur- 
faces, we note that the term 4Sr /9^/^^V arises because 
molecules within a layer of volume Vb ~ Sr can on av- 
erage reach the pore surface within time t. Following 
Ref. ^ and allowing for dg > L, it is easy to see that 
Vb cc {dg/Lf{L/r)^'r^/{i~Ds). Hence Eq. |l| becomes 
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where ^ is a constant. This expression makes it clear 
that the analysis of Ref. |l| requires L cc dg, but also r to 
be independent of dg. For a given sample with fixed dg 
and L, the limiting form of D{t) at short times is 



l-i?(i)/i?oocr3-^-' oct(3-^=)/2. 



(3) 



To illustrate this method, we show some preliminary 
PFG NMR data on Indiana limestone and a packing of 
plastic beads; experimental details are given elsewhere 
Q. Figure [| shows a log- log plot of 1 — D{t)/Do vs. r. 
For a fractal surface, the data should fall on a straight 
line with slope 3 — Dg- As seen in previous studies, the 
data for smooth plastic beads are consistent with Dg = 2. 
However, the limestone data fall on a line with distinctly 
smaller slope, giving Dg — 2.58 ± 0.14. Optical micro- 
graphs of this limestone reveal that the pore surface has 
a wide range of features on scales above 1 /im. Although 
these data span less than a decade of length scales, the 
difference in slope between limestone and plastic beads 
is unmistakable. If the grain surfaces studied in Ref. |l| 
were truly fractal, they would show D(t) time depen- 
dence similar to the limestone data in Fig. 0. 



FIG. 1: Time-dependent diffusion data for water-saturated 
samples of Indiana limestone (squares) and unconsolidated 
15 fim dia. polystyrene beads (circles). The lines show least- 



squares fits to Eq. 
at 2. For the solid 



3. For the dashed line Dg was fixed 
ines Ds was allowed to vary, yielding 
Ds = 2.58 ±0.14 for Indiana limestone and Ds = 2.2 ±0.4 for 
polystyrene beads. To accurately determine these power laws 
it was necessary to separately measure Do on a bulk water 
sample in the same apparatus, at the same temperature. 
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